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[57] ABSTRACT 

An assembly planning system that operates based on a 
recursive decomposition of assembly into subassem- 
blies, and analyzes assembly cost in terms of stability, 
directionality, and manipulability to guide the genera- 
tion of preferred assembly plans. The planning in this 
system incorporates the special processes, such as clean- 
ing, testing, labeling, etc. that must occur during the 
assembly, and handles nonreversible as well as revers- 
ible assembly tasks through backward assembly plan- 
ning. In order to increase the planning efficiency, the 
system avoids the analysis of decompositions that do 
not correspond to feasible assembly tasks. This is 
achieved by grouping and merging those parts that can 
not be decomposable at the current stage of backward 
assembly planning due to the requirement of special 
processes and the constraint of interconnection feasibil- 
ity. The invention includes methods of evaluating as- 
sembly cost in terms of the number of fixtures (or hold- 
ing devices) and reorientations required for assembly, 
through the analysis of stability, directionality, and 
manipulability. All these factors are used in defining 
cost and heuristic functions for an AO* search for an 
optimal plan. 


13 Claims, 6 Drawing Sheets 
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(DEFSCHEMA L2 
((INSTANCE LIAISON) 

(CLASS FLOATING) 

(MATING PARTS P2 P3) 
(INTERCONNECTION ATTACH) 
(FEATURES(F1.P2 F1.P3(+X -X))) 
(CONNECTION TYPE PLACE ON))) 
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(DEFSCHEMA P3 (DEFSCHEMA F2.P14 

((INSTANCE PART) ((INSTANCE FEATURE) 

(CLASS GENERAL) (FEATURE OF P3) 

(FEATURE-LIST FI .PI 4 F2.P14) (TYPE CONNECTION) 
(COORD(7 0 10 0 90 0)) (SHAPE(CYLINDER 6 1.8)) 
(VOLUME((CYLINDER 14 1) (LOCATION^ 0 0 0 0 0)))) 
0 0 0 0 0 0 ))) 
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BACKWARD ASSEMBLY PLANNING WITH DFA 
ANALYSIS 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the contractor has elected not to retain 
title. 10 

This application is a continuation of application Ser. 
No. 07/942,499, filed Aug. 28, 1992 now abandoned. 

BACKGROUND OF THE INVENTION 

L Technical Field 15 

The invention relates to a method and apparatus for 
performing backward assembly planning and operation 
using backward search of assembly sequences. 

2. Background Art 

Assembly planning can be based either on the for- 20 
ward search of assembly sequences or on the backward 
search of assembly sequences, given geometric and 
topological information on the product to be assembled. 
The backward search of assembly sequences is often 
implemented in the form of disassembly planning, as- 25 
suming that ah assembly sequence can be obtained by 
reversing the order of the corresponding disassembly 
sequence. Although there are cases where assembly 
sequences may not be the reverse of the corresponding 
disassembly sequences, disassembly planning is widely 30 
used in the research community due to its advantage in 
planning efficiency: The disassemblability of a part or a 
subassembly directly implies the satisfaction of prece- 
dence relationship, whereas, in the forward search, the 
satisfaction of precedence relationship between a pair of 35 
mating parts may not be known immediately until an 
exhaustive search is completed. 

Much effort has been expended to make assembly 
planning or disassembly planning more autonomous, 
more efficient, and closer to reality. Early work on 40 
interactive assembly planning was concerned about 
formulating a necessary and sufficient set of questions to 
be answered by the human designer that; leads to the 
complete specification of precedence relationships with 
the minimum number of question-and -answering opera- 45 
tions. Subsequently, the emphasis has been moved to 
automatic reasoning of geometric interference or path 
planning as a means of identifying precedence relation- 
ships in assembly. The feasibility of automating geomet- 
ric reasoning in assembly planning is fueled by the de- 50 
velopment of powerful computer-based geometric 
modeling and reasoning tools. The automatic test for 
the existence of a mating path in disassembly or assem- 
bly is still computationally very expensive; and this 
significantly limits the number of parts that can be han- 55 
died by an automatic assembly planner, since the num- 
ber of feasible assembly sequences may grow exponen- 
tially, so does the number of required tests for path 
existence. One of the current research goals is thus to 
achieve planning efficiency by avoiding unnecessary 60 
tests for path existence, as proposed by such methods as 
the logical inference of path existence based on the 
previous results on path planning, the test of necessary 
conditions for successful mating, which are simple to 
implement, prior to the application of path planing, the 65 
use of path finding heuristics. 

It becomes clear that, to bring assembly planning 
closer to reality, we need to consider not only identify- 
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ing the geometric interference in part mating, but also 
many other factors which affect assembly planning in a 
real factory environment: special processes such as 
testing, labeling, painting, grinding, etc. to be per- 
formed during assembly; physical constraints due to 
magnetic and electrostatic interference; stability of 
parts and subassemblies which directly affect assembly 
cost due to required fixtures; part and subassembly ma- 
nipulability; the required number of reorientations dur- 
ing assembly; tolerance propagation; as well as the 
physical constraints imposed by assembly layouts. The 
above factors should eventually be incorporated into 
the process of identifying precedence relationships as 
well as selecting preferred assembly sequences. For 
instance, the required special processes, the physical 
interference, and the stability requirement can be used 
in precedence ^identification together with the geomet- 
ric interference condition. While, the stability, the ma- 
nipulability, the assembly direction and reorientation, 
and the tolerance associated with parts, subassemblies, 
and their matings can be used in selecting prefered as- 
sembly sequences. However, the selection of a cost-ef- 
fective assembly sequence is not simple and remains as 
one of major research issues in assembly planning. This 
is not only due to the fact that evaluating the above 
mentioned conditions is difficult and computationally 
expensive but also the optimality in assembly planning 
may not be achievable through the local optimization of 
the above conditions but requires the consideration of 
global effect. 

The present invention includes a process for back- 
ward assembly planning and performance which takes 
the above issues into consideration. 

SUMMARY OF THE INVENTION 

A method of selecting and performing a sequence of 
sub-assembly operations to construct an assembly be- 
gins by identifying all direct subassemblies of the assem- 
bly by identifying potential subassemblies of the assem- 
bly, and, for each potential subassembly, determining 
whether the potential subassembly is a direct subassem- 
bly of the assembly by: (1) determining whether there is 
a common axis of separation between the potential sub- 
assembly and the remainder of the assembly, (2) deter- 
mining whether the potential subassembly is stable, and 
(3) determining whether there is a direct unobstructed 
path between the potential subassembly and the remain- 
der of the assembly. The identifying potential subassem- 
blies is performed by identifying all liaisons connecting 
pairs of nodes of each potential subassembly of the 
assembly, and, for each of the liaisons, merging liaisons 
for which there is no accessibility for interconnection 
between the subassembly and the assembly and merging 
liaisons which have no force-deliverable paths to the 
assembly. It is further performed by identifying special 
processes of the assembly of the type comprising testing 
and cleaning which require certain nodes of the assem- 
bly to be processed together during a special process. 
The identifying potential subassemblies further consists 
of grouping together nodes of the assembly which must 
be processed together during special processes of the 
type comprising cleaning and testing of subassemblies, 
whereby to prevent identification of potential subas- 
semblies comprising corresponding nodes of the group- 
ing. The overall process further consists of identifying 
precedent constraints corresponding to required se- 
quences of assembly of parts of the assembly, and the 
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identifying potential subassemblies further consists of 
grouping together nodes of the assembly which must be 
assembled in a hierarchical sequence corresponding to 
the precedent constraints, whereby to prevent identifi- 
cation of potential subassemblies comprising corre- 5 
sponding nodes of the hierarchical sequence. The deter- 
mining whether potential subassembly is stable consists 
of determining that one of (a) the potential subassembly 
has more only one translational degree of freedom and 
only one rotational degree of freedom, the one rota- 10 
tional degree of freedom being about the translational 
axis only, or (b) any floating of disconnected liaison of 
the potential subassembly includes a node of the remain- 
der of the assembly. The determining whether the po- ^ 
tential subassembly is stable consists of determining that 
the potential subassembly is unstable if either the poten- 
tial subassembly is unstable or the remainder of the 
assembly is unstable. The process further consists of 
computing a cost of mating and assemblying each direct 2 q 
subassembly, by determining the number of part reori- 
entations required to construct the subassembly and 
determining the number of mating operations required 
to construct the subassembly. The process further con- 
sists of representing the assembly as a tree of AND and 25 
OR nodes beginning with all direct subassemblies as 
bottom nodes of the tree and concluding with the as- 
sembly as the top node of the tree, the tree comprising 
branches connecting adjacent ones of the nodes, the 
branches corresponding to alternative decompositions 30 
of the assembly and computing a local cost of each OR 
node of the tree and eliminating branches correspond- 
ing to OR nodes of higher costs, and performing an 
assembly operation of subassemblies corresponding to 
the tree after the eliminating. The local cost computa- 35 
tion consists of determining the relative cost due to the 
stabilization of the direct subassembly and the remain- 
der of the assembly by using holding devices, and deter- 
mining the relative cost due to reorientation operations. 

40 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an exploded view of a flashlight assembly 
employed as a tutorial example of an application of the 
present invention. 

FIG. 2 is a liaison graph of the flashlight assembly of 45 
FIG. 1. 

FIG. 3 illustrates the special process forest of the 
flashlight assembly of FIG. 1. 

FIGS. 4 a and 4 b illustrate two assemblies in which an 
assembly sequence is not the reverse of a disassembly 
sequence. 

FIG. 5 illustrates the clustering process of the flash- 
light assembly of FIG. 1. 

FIG. 6 illustrates the special process forest of FIG. 3 55 
after removing the root node. 

FIGS, la and lb illustrate how the internal freedom 
of motion of a subassembly is determined. 

FIG. 8 illustrates an AND-OR tree employed in the 
invention as the search space for the AO* algorithm. 50 

FIG. 9 illustrates the AND-OR tree corresponding to 
the flashlight assembly of FIG. 1. 

FIG. 10 illustrates the merging and grouping opera- 
tions corresponding to the liaison graph of FIG. 2. 

FIG. 11 illustrates the overall process of the inven- 65 
tion. 

FIG. 12 illustrates a system for carrying out the in- 
vention. 


DETAILED DESCRIPTION OF THE 
INVENTION 

1 Assembly Representation 

An assembly A is a cluster of parts assembled to- 
gether by a certain assembly sequence, which maintains 
a particular geometric relationship among parts. More 
specifically, a cluster of parts, P, is considered as an 
assembly, iff 

1) Every part in P has at least one geometric con- 
straint in relation with other parts in P. 

2) There exists at least one assembly sequence which 
can put all the parts of P together to satisfy the 
geometric constraints imposed on themselves. 

3) P is stable in a sense that it can maintain its geomet- 
ric constraints either by itself or through the aid of 
holding devices. 

The above definition extends the notion of conven- 
tional assembly to include a disconnected and non-self- 
stable assembly that may be generated during the assem- 
bly of the product. 

Formally, an assembly A is represented by 4-tuples : 
A = {P(A),Gl(A),G/<A), (A)}, 

where P(A) represents the set of parts constituting A, 
Gl(A) the attributed liaison graph representation of A, 
GAA) the special process and constraint forest associ- 
ated with A, and (A) the set of all the feasible assem- 
bly sequences for A. 

1.1 Attributed Liaison Graph 

A liaison is said to exist between a pair of parts if one 
part constrains the freedom of motion of the other by a 
direct contact. A liaison graph is a graph, G, G=(N, E), 
with N representing a set of nodes, and E representing 
a set of edges. A node n, n N, is assigned to each part 
of an assembly, and an edge e, e E, is assigned to a 
liaison between a pair of parts. An attributed liaison 
graph is a liaison graph with frames attached to individ- 
ual nodes and edges of the liaison graph to describe the 
attribute associated with a node or an edge. A part 
frame attached to each node describes the attributes 
associated with a part, including 1) the part geometry 
and coordinate, 2) the mating volumes and the contact 
surfaces as part features, and 3) the physical properties 
of a part such as weight. An edge frame attached to 
each liaison describes attributes associated with the 
liaison. The attributes of a liaison consist of 1) the parts 
in associated with the liaison, 2) the mating features and 
mating directions, 3) the mating type such as insertion 
and place-on, and 4) the interconnection type such as 
Attachment, Force-fit, Connectors/Retainers, Push- 
and-Twist, Screw, Glue, or Welding. Thus, Gl(A) 
contains information on the topology of part configura- 
tions, the geometry and relative pose of parts in A, the 
interconnection mechanisms of part connections, and 
the local freedom of motion in part mating. 

1.2 Special Process and Constraint Forest 

The assembly of a product involves not only the 
interconnection of parts to form required liaisons but 
also the execution of special processes such as testing, 
cleaning, surface treatment, painting and packaging 
during assembly. In addition, there may exist assembly 
constraints that should be satisfied in order to prevent 
electrostatic, electromagnetic, and thermal interfer- 
ence, as well as mechanical vibrations and chemical 
pollutions, during assembly. The accomplishment of 
special processes and the satisfaction of assembly con- 
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straints impose precedence relationships in assembly. 
Especially, assembly precedence relationships should 
be generated in such a way as to guarantee the execu- 
tion of multiple special processes in a certain prece- 
dence order. 5 

Special processes and assembly constraints associated 
with A are represented by the special process and con- 
straint forest, G^(A), consisting of the special process 
forest, Gs(A), and the assembly constraint forest, G C (A). 

The special process forest, G*(A), consists of a set of 10 
trees having the following properties: 

1) A node n^of G^A) represents a special process, S/. 
n/S is associated with a tuple (Pf, 2/*) where P/* 
represents a set of parts involved in S /, and 2/* 
represents the union of the parts involved in S/and 15 
all the special processes that should precede S/: 
2/=P/U{Uy2/}, V n/: n/is a child of n/*(a recur- 
sive definition). 

2) A branch b/ 5 connecting the two nodes, n if and n/, 
represents the precedence relationship between S/ 20 
and Sf. Sy<S/if n/is a child of n/. Special processes 
corresponding to the sibling nodes of a tree or to 
the nodes of different trees have no precedence 
relationship. 

The assembly constraint forest, Gc(A), consists of a 25 
set of trees having the following properties: 

1) A node n^of Gc(A) represents a liaison 1/of Gi(A). 
n,* is associated with a tuple (P^, 2^), where P,* 
represents a pair of parts involved in 1/, 2,* repre- 
sents the union of the parts involved in 1/and all the 30 
liaisons that should precede If. 2/ c =P/ c U{Uy 2/}, 

V n/, n f is a child of n^. 

2) A branch bf connecting nf and n f represents a 
precedence relationship between 1/and 1/. ly<l/if nf 

is a child of n/°. Liaisons corresponding to the sib- 35 
ling nodes of a tree or to the nodes of different trees 
have no precedence relationship. 

Example: Flashlight Assembly 

The flashlight assembly shown in FIG. 1 consists of a 
body(BD), two dry cells(Bl,B2), a bulb carrier(BC), a 40 
bulb(BL), a receptacle(R), lens(G), and a cover(C). 
FIG. 2 shows the attributed liaison graph representation 
of the flashlight assembly. FIG. 3 illustrates the special 
process forest associated with the flashlight assembly. It 
shows that the cleaning of the receptacle(R) should 45 
precede the testing of the bulb unit(BC,BL), such that 
BC, BL, and R should be clustered in a subassembly 
during backward assembly planning. 

2 Backward Assembly Planning 

2. 1 Definition 50 

A subassembly of an assembly A, S/| A, is an assembly 
that can be generated in one of the feasible assembly 
sequences for A. A direct subassembly of A, S/4| A, is 
defined as a subassembly of A, which can be directly 
assembled with A— S/4| A at the last step of one of the 55 
feasible assembly sequences for A. The assembly plan- 
ning of an assembly A, A={P(A), G^A), Gp(A) f 
11(A)}, is the process of generating a set of assembly 
sequences Ilc(A)Cn( A ) based on the given P(A), 
Gl(A) and Gp(A), and the criteria for selecting pre- 60 
ferred assembly sequences. The backward assembly 
planning of A, BAP(A), represents a particular method 
for achieving the assembly planning of A, based on the 
recursive identification and selection of preferred direct 


6 

i= 1, . . . , m}). Then, it recursively applies the process 
of decomposition to the subassemblies generated by the 
previous decomposition until no further decomposition 
can be applied. 

2.2 Backward Assembly Planning vs. Disassembly Plan- 
ning 

BAP(A) differs from disassembly planning in that 
BAP(A) can handle the case where an assembly se- 
quence can not be obtained from the reverse of a disas- 
sembly-sequence. For instance, a sequence of opera- 
tions to disconnect a liaison of force-fit interconnection 
type is often quite different from the reverse of a se- 
quence of operations to interconnect the liaison. As 
shown in FIG. 4a, the snaps of part A should be wid- 
ened in order to disassemble part B from part A, which 
may require a new tool. A more interesting example is 
given in FIG. 4 b, where we can have the following 
disassembly sequence: C-A-D-B, since Screw C can be 
disassembled before Cover A is removed. However, the 
reverse of the above disassembly sequence, B-D-A-C, 
can not be an assembly sequence, since the inability of 
holding part D during assembly of screw C prohibits 
such an assembly sequence. BAP(A) is able to identify 
that the direct subassembly of (A,B,C,D) is A but not C, 
and generate a feasible backward assembly sequence, 
A-C-D-B. 

2.3 Identification of S/*| A 

The core of BAP(A) is the identification of Sf*\ A. 
For the identification of S/4| A, we introduce the fol- 
lowing conditions that a cluster of parts P | A, P | A'C A, 
can be a direct subassembly of A: 

1) Accessibility Condition: P|A is accessible by a 
gripper or a tool. 

2) Stability Condition: P|A and A— P| A are stable, 
i.e., they either are self-stable or can be stabilized 
by external devices. 

3) Local Mating Motion Condition: All liaisons be- 
tween P | A and A — P | A have at least one common 
axis of separation. 

4) Path Existence Condition: P | A can be brought to 
A — P| A from the free space for melting. 

5) Interconnection Feasibility Condition: P | A can be 
interconnected to A— P|A by applying the inter- 
connection operations defined for the liaisons be- 
tween P| A and A — P| A. 

6) Process Constraint Condition: P | A meets the con- 
straints defined by the special process and con- 
straint forest of A, Gp( A). 

To avoid the unnecessary evaluation of computation- 
ally expensive conditions, such as the path existence 
condition, we hierarchically evaluate the above condi- 
tions. The following sections describe the details of the 
method for testing the above conditions. 

3 Part Clustering 

The identification of S/4 [A can be accomplished by 
checking the above conditions for individual P|A’s 
defined by the cut-sets of G/,(A). Since the number of 
cut-sets of Gl(A) is often large and grows exponentially 
with respect to the number of nodes in G//A), testing 
all the cut-sets of G^(A) results in inefficiency in assem- 
bly planning (especially due to the high cost of testing 
the path existence condition). 

The number of the cut-sets that require the test of 
path existence condition may be reduced considerably if 


subassemblies. BAP(A) first identifies and selects a di- 65 we first select those cut-sets that satisfy the interconnec- 
rect subassembly of A, S/|A(or a set of direct subas- tion feasibility condition, the process constraint condi- 
semblies of A, {S/4| A, i= 1, . . . , m}) and decomposes tion, and the stability condition, prior to the evaluation 

A into S/4|A and A — S,4| A (or {S^| A and A— S/4] A, of the path existence condition. For the selection of 
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those cut-sets satisfy the interconnection feasibility and 
process constraint conditions, we first construct an 
abstract liaison graph, Gl(A), from Gi,(A), as follows: 

1) Merging those parts of A that cannot be intercon- 
nected, although it is assumed that those parts are 
already brought into their mating position. This 
process identifies the liaisons of G L(A) that violate 
the interconnection feasibility condition, by rea- 
soning whether the preconditions for a particular 
interconnection mechanism associated with a liai- 
son can be satisfied. 

2) Grouping those parts that should belong to the 
same subassembly to meet the requirement of spe- 
cial processes. 

The parts merged together become a supernode, 
whereas the parts grouped together become a group 
node in Gl(A). Then, the cut-sets of Gl( A), where a 
super node and a group node are considered as a single 
node, satisfy the interconnection feasibility and process 
constraint conditions. 

3.1 Interconnection Feasibility and Part Merging 

A liaison 1/ between P|A and A— P|A cannot be 
completed, even if P | A can be brought into its mating 
position, if the following conditions are not satisfied: 

1) The feasibility of applying tools and installing con- 
nectors required for the interconnection of 1 /. 

2) The feasibility of applying a force, while maintain- 
ing stability, required for the interconnection of 1 /. 

The test of the first condition can be done by check- 
ing the existence of an open channel to the designated 
part locations, through which tools and connectors can 
be operated without geometric interference. 

The test of the second condition requires reasoning 
on the force delivery to a liaison 1 / through intermediate 
liaisons. To be more specific, let us first categorize a 
liaison into one of the following three classes: 

1) A liaison is said to be floating if there exists no 
physical force holding the parts (associated with 
the liaison) together. For instance, a liaison with 
the interconnection type of “attachment” is a float- 
ing liaison. A floating liaison may or may not be 
stable, depending on the geometric constraints 
imposed on the parts associated with the liaison. 

2) A liaison is said to be rigid, if there exists physical 
force holding the associated parts together, by 
which the liaison becomes self-stable even under 
the presence of external force. For instance, a liai- 
son with the interconnection type of “force-fit”, 
“welding”, or “connectors” may be classified as a 
rigid liaison. 

3) A liaison is said to be firm if there exists a physical 
force holding the associated parts together, by 
which the liaison becomes self-stable, although it 
may exhibit a deformation or a freedom of motion 
under the presence of external force. 

A liaison is associated with its local freedom of mo- 
tion: 

The local freedom of motion of a liaison 1/, LFM(h), 
connecting the two parts, Pi and P 2 , is represented by 
the freedom of motion of Pi against P 2 , LFM(lf, Pi | P 2 ) 
or the freedom of motion of P 2 against Pi, LFM(1/; 
P 2 |Pi), where LFM(1,; Pi|P 2 ) or LFM( 1 , ; P 2 |Pi) is 
represented in terms of the coordinates of the assembly 
to which Pi and P 2 belong. 

A local freedom of motion of a liaison may be 
changed after the interconnection is completed: A rigid 
liaison completely loses its local freedom of motion 
after the interconnection is completed. A firm liaison 
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may show a local freedom of motion, when the force 
applied to the liaison is more than a certain threshold. 
To distinguish the local freedom of motion of 1/ after the 
interconnection from that before the interconnection, 
LFM + (1/) is used to represent LFM(1/) after the inter- 
connection. Note that, due to the orthogonality be- 
tween motion space and force space, a liaison connot 
deliver force to the direction where a local freedom of 
motion exists. 

To test the feasibility of applying force required for 
the interconnection of 1 /, let us define the following: 
Definition: An Accessible Node, A-node 

A part is accessible by a tool if it is reachable and 
graspable by a tool for an assembly operation. A node is 
accessible if any of the parts forming the node is accessi- 
ble. 

Definition: An Access Path to a node n, A-path 

An access path to a node n, A-path, is a path starting 
from an A-node and ending with the node n without 
having any other A-node in the middle of the path. By 
definition, an A-node has an A-path to itself. An access 
path is represented by an ordered set of liaisons or parts 
on the patch. A node may have one or more A-paths. A 
pair of A-paths, A-path(n^) and A-path(n 2 ), n 1 = 7 ^ 12 , are 
said to be independent each other, if they share no com- 
mon part. 

Definition: The Internal Motion Space, M [A-path(n)], 
and the Static Force Space, T [A-path(n)], associated 
with A-path(n) 

The internal motion space, M [A-path(n)], of A- 
path(n) is defined by the union of the local freedom of 
motion of individual liaisons in A-path(n), and repre- 
sents the flexibility that the configuration of parts along 
A-path(n) can be deformed by an external force, with 
the first part (corresponding to A-node) and the last 
part (corresponding to the node n) fixed in space. Note 
that M, [A-path(n)] is a function of the external force 
applied to A-path(n) when A-path(n) includes firm liai- 
sons, since the external force determines which firm 
liaisons can be broken. As will be explained shortly, the 
external force to be applied to A-path(n) is given as the 
force required for the interconnection of the liaison 
which n is associated with, and is subject to the evalua- 
tion of force-deliverability. 

Assume that A-path(n) is represented by an ordered 
set of liaisons, {li, I 2 , . . . , l r }, with 1 / formed by a pair of 
nodes, (n/i, n/ 2 ), and that (n/j, n/ 2 ) is ordered along A- 
path(n) in the direction toward n. Then, M [A-path(n)]- 
= U /= i r LFM(l z ; n/i | n/ 2 ). The static force space, JT[A- 
path(n)], of A-path(n) defines the static force that can be 
delivered to the node n through A-path(n). J[A- 
path(n)] is represented by the orthogonal complement 
of M [A-path(n)], i.e., ^[A-path(n)] = {±x, ±y, ±z, 
±vff, ±0, U/=/LFM (lr, n/i|n/ 2 ). 

Definition: A Force-Deliverable A-path 

A-path(n/i) is said to be force-deliverable to n/i for 
the liaison 1/ of (n/i, n/ 2 ) if F[A-path(n;i)] includes the 
force required for the interconnection of n/j to n/ 2 . 

The interconnection feasibility condition of a liaison 
of (n/i, n/ 2 ) is now transformed to the verification of the 
existence of an independent force-deliverable A-path 
for n/i and n/ 2 . The force-deliverability of an A-path to 
n/i or n /2 depends on the force required for the intercon- 
nection of 1/. That is, M [A-path(n/i)] or M [ A-path(n/ 2 )] 
is determined by the freedom of motion of individual 
liaisons along the A-path under the presence of an exter- 
nal force equivalent to the force required for the inter- 
connection of 1/. For instance, in case 1/ is floating, the 
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amount of force required by 1/ is negligible, and the 
decision of the internal motion space of the A-path is 
based solely on the floating liaisons along the path. Note 
that, in this case, the force-deliverability of 1/ becomes 
equivalent to the feasibility of maintaining stability dur- 5 
ing interconnection. 

As a summary, we present the following merging 
principle: 

Merging Principle: 

The liaison 1/of (n/i, n/ 2 ) can not be listed as a cut-set, 10 
and consequently n/j and n/ 2 should be merged together, 
if one of the following conditions is true: 

1) It is not feasible for the tools and connectors (re- 

quired for the interconnection of 1,) to access the 
designated locations. 15 

2) Either n/j or n/ 2 has no independent force-delivera- 
ble A-path, including the case where either n/i or 
n/2 has no A-path at all. 

The following algorithm summarizes the process for 
part merging: 20 

Algorithm: Part Merging 

Input: Gl(A). 

Output: Gl(A) with super nodes. 

Method: 

Step 1. Put all the liaisons of Gl(A) in Open set. 25 
Identify A-nodes of Gz,(A ). 

Step 2. If Open set is empty, then stop. 

Step 3. Select a liaison 1/ from Open set in an order 
starting from an A-node, and remove L/from Open 
set. 30 

Step 4. Check whether 1/ requires tools or connectors 
to complete the interconnection. If not, go to Step 
6 . 

Step 5. Check the accessibility of tools or connectors 
to the designated locations. If not, merge n/i and 35 
n/ 2 associated with 1/ and go to Step 2. 

Step 6. Check whether n/j and n/ 2 have independent, 
force-deliverab A-paths. If not, merge n*/i and n/ 2 
and go to Step 2. , 

Example: Part Merging of the Flashlight Assembly 40 

The A-nodes of the flashlight assembly are identified 
as BD and C (refer to FIG. 2). To show how each 
liaison in FIG. 2 can be tested for the merging, let us 
consider I9 as an example: 

For I9 of (C,R), A-path for C is {C}, whereas A-paths 45 
for R are {C,R}, {C,G,R>, {BD,B2,BC,R}, 

{BD,B1,B2,BC,R} and {BD,B1,B2,BC,BL,R}. 1 9 does 
not require a tool or a connector to complete the inter- 
connection. The part, C, has a force-deliverable A-path, 
{C}; hpwever the part, R, does not have a force-deliv- 50 
erable A-path independent of {C}. This is because all 
A-paths for R from BD can not deliver force to R due 
to the floating liaisons li, I2, U, I5, and I7. Therefore, I9 is 
merged. 

FIG. 5 shows the result of the merging process. In the 55 
drawing, a circle represents a node; a dash line repre- 
sents a floating liaison; a thin solid line represents a firm 
liaison; a heavy solid line represents a rigid liaison; sin- 
gle hatching represents merging; and cross-hatching 
represents grouping. 60 

3.2 Special Process Constraints and Part Grouping 

The recursive determination of S/^|A in backward 
assembly planning should support the execution of spe- 
cial processes in an order specified by G^A) and satisfy 
the assembly constraints specified by Gc(A). This re- 65 
quires the identification of those parts of A that should 
be grouped together, or should not be divided into 
S z ^| A and A — S,^| A, in order to implement the special 
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process forest G^A) and the assembly constraint forest 
G C (A) as follows: 

1) Given G^A), we first determine which special 
processes should be executed with A, i.e., the selec- 
tion of one or more special processes executable 
with A from G/A) in a top-down order from the 
root nodes (note that there can be multiple trees in 
G^A) ). We then remove out those nodes corre- 
sponding to the selected special processes from 
G^(A), resulting in G^A). Since the special pro- 
cesses remaining in G'/A) need to be accomplished 
in the later stage of backward assembly planning, 
we should group those parts in the accumulated 
part list of each of the root nodes of G'^A). 

2) In order to preserve the precedence relationship 
defined by G C (A), we should group those parts in 
the accumulated part list of each child of individual 
root nodes of G C (A). Note that a liaison repre- 
sented by a root node of G C (A) is eligible for de- 
composition at the current stage of backward as- 
sembly planning, and is exempt fibre grouping. 

Grouping Principle: 

For a given A, G/,(A), Gs(A) and G C (A), we group 
those nodes of G^(A) that belong to either of the fol- 
lowing lists: 

1) 2o i s associated with each root node noA of G'XA) 
where G'^A) is obtained by removing out from 
G*(A) those nodes that are selected for processing 
with A. 

2) disassociated with each child, n\f, of individual 
root nodes of Gc(A). 

Example: Part Grouping of the Flashlight Assembly 

First, G'XA) is obtained in FIG. 6, since the testing of 
the whole assembly should be done with A. 

Then, {B1,B2} and {BC,BL,R}, the accumulated 
part lists of individual root nodes of G's(A), become 
groups and .{BC,BL,R} becomes another group, as 
shown in FIG. 5. {G} is a single part and is not indeed 
in the group list. 

4 Identification of Direct Subassemblies 

The problem of finding direct subassemblies is now 
transformed into the problem of finding valid cut-sets of 
the abstract liaison graph. Gx(A). A cut-set, y/, decom- 
poses Gx(A) into disjoint subgraph, Gl Si (A, y /) and 
G l 52 (A, y/) where Gx 5l (A, y /) or Gx^A, y /) may not 
be a connected graph, but may be a collection of multi- 
ple connected subgraphs, Gx 5l (A, y/)=={6i:‘ s V(A, y/), 
j=l, 2, . . . , li)} or Gr*2(A, y/)={g^y(A, y/), j = l, 2, 
. . . , 1 2 }. A cut-set, y /, is valid if: 1) Each connected 

subgraph, Gx si -/(A, y /), j = l, 2 li, Gx S2 -/(A, yi), 

j= 1, 2 . . . , 1 2 , generated by the cut-set, y/ includes at 
least one A-node. 2) There exists a collision free path 
between the cluster of parts, P|A corresponding to 
Gl S i ( A, yi) and A— P | A, corresponding to Gx^A, yi) 
for their mating. 

The first condition ensures that P|A and A— P|A, 
can be handled by a tool for assembly, whereas the 
second condition ensures that S/*| A can be placed in its 
mating position. For the test of the second condition, let 
us define the following: 

1) The predicate LM(P| A, A— P] A)is true if all the 
liaisons between P | A and A — P | A has at least one 
common axis of separation. 

2) The predicate PE(P| A, A— P| A) is true if there 
exists a path along which the cluster of parts P | A 
can be brought to its mating position without col- 
liding with the rest of the assembly. 
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Since the premise that LM(P | A, A— P | A) is not true 
is a sufficient condition that PE(P|A, A— P| A) is not 
true, the test of LM(P|A, A — P|A) preceding the test 
of PE(P | A, A— P | A)can provide a considerable reduc- 
tion in the number of the costly PE(P | A, A — P | A) test. 5 
Algorithm^ Identification of Direct Subassemblies 
a Input: Gl(A), A-Node Set {all the A-nodes of 
Gl(A)}, Tested Cut-set List={0}, Valid Cut-set 
List = {0}. 

Output: A list of direct subassemblies specified in 10 
Valid Cut-set List. 

Method: 

Step 1. If A-node Set is empty, stop. 

Step 2. Select an A-node from A-node Set, and re- 
move it from A-node Set. 15 

Step 3. If there exists^ cut-set %, y/, y /^Tested Cut- 
set List, such that Gl Si £A, y /) includes the A-node 
selected in Step 2 and Gr£ 52 (A, y ; ) includes at least 
one A-node of Gl(A), then continue. Otherwise, 
go to Step 1. 20 

Step 4. Put y/into Tested Cut-set List. 

Step 5. Test y/for LM. If it is false, go to Step 3. 

Step 6. Test y/for stability. That is, if one of the subas- 
semblies produced by y/ is unstable or the subas- 
semblies produced by y/ have no common assembly 25 
pose among themselves go to Step 3. 

Step 7. Test y / for PE. If it is true, put y,- in Valid 
Cut-set List. Otherwise, go to Step 3. 

The details of the method for testing the stability of 
y,- will be presented in the next section. 30 

5 Subassembly Evaluation 

The direct subassemblies identified from Gl(A) are 
subject to further evaluation for the selection of pre- 
ferred direct subassemblies. 

The evaluation of a subassembly is based on the fol- 35 
lowing criteria: 1) Stability, 2) Directionality, and 3) 
Manipulability, which leads to an estimate of assembly 
cost involved in local assembly through the identifica- 
tion of the required holding devices and reorientations 
during assembly. Note that the optimal selection of a 40 
direct assembly at each stage of backward assembly 
planning based on local evaluation may not yield a 
globally optimal plan. As will be shown later, this prob- 
lem is handled by the AO* algorithm with its cost and 
heuristic functions defined in terms of the above crite- 45 
ria. 

5.1 Stability 

The stability of a subassembly can be represented in 
terms of the internal freedom of motion associated with 
the subassembly. To present a method formalizing the 50 
internal freedom of motion of a subassembly, let us first 
define the following: A cluster of parts in S^\A, 

Pa | (S^| A). is said to be floating if it ^connected to the 
rest of Sf* | A only by floating liaisons. P a | (S t d | A) corre- 
sponds to a subgraph of Gl(S,^|A) that can be sepa- 55 
rated from Gl(S^| A) by a cut-set consisting only of 
floating liaisons (denoted here as a floating cut-set). A 
cluster of parts in S/|A, Pa|(S^| A), is said to be dis- 
connected, if it has no liaison connected to the rest of 
S/rf | A. Pa ] (S/^| A) corresponds to a disjoint subgraph in 60 
Gz,(S^|A). The stability of S^| A can be defined based 
on a set of floating clusters of # parts, Pa|(S/^|A), and 
disconnected clusters of parts, Pa|(S, j |A), included in 
S,^|A. Let us assume that a floating cut-set, r*, of 
Gz,(S/f | A) decomposes S Z ^|A into P*|(S^|A) and S t d ~ 65 
| A— Pa| (S t d \ A). The local freedom of motion of Fa, 
LFM(Fa), can be defined as the local freedom of motion 
of Pa{S,^|A) against S t d \ A-P*|(S^| A): 


LFMfid £ LFM(P k \(S d \A),Sfi\A - P k \{S d \A)) 

= nLFMVj^iS^A^A - 

for all ly, IjEjk- 

Note that Pa|(S^| A) should be chosen in such a way 
that S/*| A— P&|(S/^| A) includes one of the A-nodes of 
the assembly A that will be used for holding S^|A 
during the mating of S/^| A with A — S/| A. Now let us 
define the following: 

Definition: Internal Freedom of Motion of S^\ A 
The internal freedom of motion of S Z ^|A, IFM- 
(Si d \ A), is defined as a collection of assembly directions 
to which S/*| A can be broken apart. IFM(S/^| A) can be 
calculated by the following rules: 

1) J Pa | (S/^ | A) and Pa|(S/^| A)=>IFM(S Z ^| A)=0; 

2) J PaKS^|A) but BPa|(S^|A)=>IFM(S^ 

| A)= U LFM(Fa), V Fa; 

3) 3P*|(S^|A)=>IFM(S^|A) = {±x, ±y, ±z, ±*, 

±0, ±4>}. 

As an example, let us consider a simple 2-D assembly 
shown in FIGS, la and lb . 

Since li, I2 and I3 are floating liaisons, we have that 
t r i = {ll, h}, t2~{li, I3}, rs={l2, b}. 

Assuming that S/| A is oriented with reference to the 
assembly pose of A as shown in FIG. 7, and that Pi, an 
A-node of A, is selected for grasping of S t d \ A, we have 

LFM(Fi)=LFM(P 2 UP 3 , Pi)={+X, +z}, 

LFMfo) — LFM(P 2 , PiUP 3 ) = {+z}, 

LFM(7 3 )=LFM(P 3 , PiUP 2 ) = {+x, +z>. 


Therefore, 

IFM(S/*|A)={+x, +z}. 

Based on the definition of IFM(S/^| A), we can estab- 
lish the stability condition for S/| A, as follows: 

1) S/*|A is said to be self-stable, if IFM(S/^| A) is null, 
or IFM(S/^| A) has at most a single translational 
freedom of motion with a rotational freedom of 
motion only about the axis of translation (i.e. 
Pa I (Si d | A) and S t d \A Pa|(S^|A) form a single 
peg-and-hole type mating relationship). 

2) S/*|A is said to be stable with the assistance of 
holding devices, if each Pa|(S/| A) or Pa|(S^| A) 
contains an A-node of the assembly A. This implies 
that the mating operation of S t d \ A can be stabilized 
and completed with the assistance of external de- 
vices holding and Pa|(S/^| A) of more than a peg- 
and-hole type of motion freedom against S Z ^|A 

P fc |(S<rf|A). 

3) Otherwise, S/| A is said to be unstable. A stable 
S t d 1 A, whether it requires a holding device or not, 
has one or more stable assembly poses, where an 
assembly pose is represented by assembly coordi- 
nate aligned with the direction of part stacking 
against gravity. For instance, S/^|A with null 
IFM A (denoting IFM after the incorporation of 
necessary holding devices) can have assembly pose 
of ±x, ±y and ±z. 

Si d \A with IFM^of {+x, can have an assembly 

pose of -j-x, requiring a reorientation of S t d \A to align 
+x with the stacking direction (against gravity). 
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Let us now consider the stability associated with an 
assembly operation; 

Definition: Stable Assembly Operation 

The assembly operation between S Z ^|A and A — 
S t d | A is said to be stable, if S/* | A and A — S t d | A have at 
least one common stable assembly pose defining a stable 
assembly direction. 

Now, the evaluation of S^| A in terms of stability is 
based on the stability of S/*| A and A— S^( A, as well as 
the stability of the assembly operation between S/|A 
and A— S t d \A: 

1 ) If either S[ d \ A is unstable or A— S/*| A is unstable, 
Si d \A can not be selected for a direct subassembly 
of A. 

2 ) When S/^j A and A— S/*| A have no common sta- 
ble assembly pose, S/| A can not be selected for a 
direct subassembly of A. 

3 ) Otherwise we evaluate the assembly cost incurred 
by the need to stabilize S^|A and A— S^|A as 
well as the assembly operation between S^[ A and 
A-S^|A. 

The assembly cost is directly related to the number of 
holding devices required for stabilizing S^| A and A — 
S/^| A, and the necessity of reorienting S 2 ^| A and A — 
S/j A for a stable assembly operation. The latter will be 
analyzed in more detail in the next section in relation to 
the directionality in assembly and the determination of 
best assembly poses. 

5.2 Directionality and Assembly Pose 

The directionality in assembly is another important 
factor affecting assembly cost. Locally, a stacking oper- 
ation is considered more cost-effective than a non-stack- 
ing operation. Globally, a single direction of assembly is 
preferred to multiple directions of assembly. Therefore, 
the evaluation of directionality in assembly should be 
based on both the local assemble direction between 
S^j A and A — S^]A, and the uniformity of assembly 
directions embedded in S^| A and A — S z ^| A. 

It should be noted that whether or not the local as- 
sembly direction between S z ^| A and A— S z ^| A can be a 
stacking direction depends on the choice of the mating 
pose (as one of the stable assembly poses common to 
S^| A and A— S t d \ A). However, the selection of a mat- 
ing pose between S t d \A and A— S t d \ A based solely on 
implementing a stacking operation may incur the need 
to reorient the assembly of S z ^| A and A— S^| A. The 
assembly pose of S t d \ A or A — S^JA should be chosen 
from the set of stable assembly poses of S z ^| A or A — 
S/^l A, which may differ from the selected mating pose 
between S^| A and A— S t d \ A. This also incurs the need 
to reorient S/ 1 A or A — | A, that S/^ | A or A — S i d | A 

is brought into its pose. This implies that the determina- 
tion of an assembly pose and an assembly direction 
should consider the trade-off between maximizing the 
directionality in assembly and minimizing the reorienta- 
tion of assembly pose. In principle, it is desirable to 
avoid a costly reorientation, unless the reorientation is 
required to allow many local stacking operations in the 
subsequent backward assembly planning, thus justifying 
the cost of reorientation. 

<Cost involved in the Local Mating Operation be- 
tween S t d \A and A — S Z ^|A> 

Let us first introduce the following notational con- 
ventions: 

1 ) ti, t2 ^ an assembly pose of S^| A and A— S^| A 
represented with reference to the previously deter- 
mined assembly pose of A, t*. 


2 ) {ti 5 }, {ti s } £ a set of stable assembly poses for 
S/*|Aand A-S/*|A. 

3 ) {ti2 5 } ^ a set of stable assembly poses common to 
S/*| A and A — S/*| A, i.e., {ti 2 s }={ti s }n{t 2 1 }. 

5 We can associate each pair (t^, t 2 5 ) where ti 5 E{ti 5 } 
and t 2 5 E{t 2 5 }, with the relative assembly cost, L, in- 
volved in a local mating operation. L can be determined 
based on the required number of reorientations and the 
directionality of mating operations (whether it is a 
1° stacking operation or a non-stacking operation ), as well 
as the difficulty of handling the related subassemblies. 

The reorientation of the assembly poses of S t d \ A and 
A— S/^| A during assembly becomes necessary due to: 

1 ) The need to transform ti 5 and/or t 2 5 into a mating 

15 pose, ti 2 s , ti2* £ t ^ ie case where ti^tj^ or 

t 2 5 = 7 ^tl 2 ^ 

2 ) The need to transform the selected mating pose, 
t\ 2 S i into the assembly pose of A, t* in the case 
where tn^t*. 

20 Table 1 shows the reorientations required for the 
mating between S^| A and (A— S^| A) under various 
conditions on ti^and t 2 * The directionality of the mat- 
ing operation can be tested by transforming the mating 
directions (between S t d \ A and A — S z ^| A) in terms of t* 
25 into the mating directions in terms of tj 2 J , where the 
mating directions in terms of t* are identified during the 
verification of the path existence. 

The relative assembly cost, L, involved in a local 
mating operation can now be calculated for individual 
30 (tv 5 , t 2 -0, by 

L^iyiREWft^REOin^+yiREOitn*)- 

]ao+74^0> 

35 TABLE 1 

The Reorientations required for the Local Mating Operation 
between Si^|A and A — S^lA 


The Required 

Conditions Reorientations 



ti s =t 2 s 

tl*£{tl2 s } 

l(ti2 5 —t*) 
3 (ti 5 — > 4 12 s , t 2 , ~* 

tl 2 s , ti2 s ->t*) 


tlW 

tj* €{ti2 s } or 
t 2 J €{ti2 5 } 
ti s £{ti2 s } and 
t2 5 £{tl2 5 > 

2(t ] ^ 1 2 s or t 2 5 — > 
tl2 J , ti2*-»t*) 

3 (t] s — >ti2 s , t 2 s — *> 

tl2* ti2 s -*t*) 


ti 5 =t 2 s 

trV=t*, ti*e{ti2*} 
tlVt*, 

0 

Htl2 5 -t*) 
2(ti 5 — t*, t 2 5 ->t*) 


t! w 

ti*=t* or t 2 *=t* 
ti^t* and t 2 'V : t* 

l(t 2 *— t* or tv*--***) 
2(ti*— t*, t 2 *— t*) 


50 

where REO(ti 5 ), REO(t 2 s ), and REO(ti 2 *) are binary 
functions of either 1 (when the reorientation of the 
corresponding assembly pose is required) or 0; ao and 
fio, represent respectively the normal relative assembly 
55 cost due to a reorientation and a mating motion, 71, y 2 , 
73 and 74 represent the effect of part manipulability on 
the relative assembly cost (refer to the next section for 
more detail). The best assembly poses of S,^|A, 
t*(S/|A), and A— S^JA, t*(A— S/*|A) can be deter- 
60 mined based on achieving the minimum relative assem- 
bly cost, L, due to the local mating operation. 

<Cost involved in the Assembly of S/| A and A — 
S d \A> 

Let us now consider the uniformity of assembly di- 
65 rection embedded in S t d \ A and A — | A, as a means of 

globally estimating the relative assembly cost, R, in- 
volved in the assembly of S^|A and A — S z ^| A. Since 
the exact evaluation of R can only be obtained after a 
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30 
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complete assembly plan is formulated, we indirectly 
estimate R based on the following two major factors 
contributing to R: 

1) The estimated relative assembly costs, R 0 (S^| A) 
and R 0 (A— S z ^| A), due to the number of reorienta- 
tions involved in the assembly of S^| A and A — 
S/*|A. 

2) The estimated relative assembly cost, Rj(S z ^|A) 
and R^A— S/*| A), due to the number of stacking 
and non-stacking operations involved in the assem- 
bly of S/*j A and A— S z ^| A. 

The estimation of R tf (S^|A) is based on the direc- 
tional uniformity of S/| A under the special consider- 
ation of the necessary reorientation(flip-over) of the 
bases of S/*| A, as follows: 

S/^|A is said to have m degrees of directionality in 
£(§=x, y, or z), denoted as D^(S ; ^| A)=m, in the case 
where there exist m liaisons, 1/, i = 1 . . . m, in Gl(S z ^| A) 
such that +§ or — §€{LFM(1/)} for i=l . . . m. Note 20 
that in defining the directionality of S/*|A, x, y, and z 
are referenced in terms of the assembly pose of A, and 
LFM(1 Z ) can be computed either by LFM(1 Z -, Pi, P 2 ) or 
LFM(1 Z ;P 2 , Pi), since LFM(1 Z ; Pi, P 2 )= -LFMfly; P 2 , 

Pi). 

| A is said to have the directional uniformity of r in 

denoted U^(S^| A)=r, in the case where D^(S;^[ A)/- 
Card {I* l/€Gz,(S/^[ A)}=r, where r^l. S z ^| A has the 
maximum directional uniformity in £, if S^| A has the 
maximum directionality in 

Now, let us select U^(S^| A) from A), £=x, 

y, or z} in a decreasing order until the accumulation of 
the selected U§(S Z ^| A)’s becomes greater than or equal 
to unity. Define {|*(S^|A)} as a set containing £’s 
which are associated with the selected U|(Si^|A) s. 
Then. {£*(S^| A)} represents the minimum number of 
distinctive assembly directions involved in the assembly 
of S Z ^|A, such that Card{|*(S^| A)}— 1 provides an 
estimate of the lower bound of the number of reorienta- 
tions required for the assembly of S z ^| A. A more accu- 
rate estimate of the required number of reorientations 
can be obtained by identifying the base flip-overs that 
become necessary when parts are connected to the both 
side of a base such that {LFM(1 Z - Parti, Base), VPart/ 45 
linked to base} includes both + § and — §. To be more 
specific on this point, let us define the following: A base 
node, njj, of S/*| A, n^S/^l A), is a node of Gz,(S z ^| A) 
which has the degree far greater than the average de- 
gree of a node of Gl(S z ^|A), i.e., the degree of ng(S z ^- 50 
|A)^k. Average Degree of Gl(S z ^|A), where k is a 
construct such that k>>l. A base node n#(S z ^|A) is 
said to have the directionality, D«fl(S/*|A) where 
D«fl(S/*|A) is the set of independent directions in- 
volved in {LFM(lf, Pi, Base), V 1/: 1/ is linked to n^}, 55 
where D„b(S^\ A) C{+x, -x, +y, — y, + 2 , -z}. In 
the case where D„^(S Z ^| A) contains -f £ and — £ (|=x, 
y, or z), it is necessary to flip-over the base to accomo- 
date + § and — 

Considering the above R 0 (S Z ^| A) can be estimated by 
the following equation: 


35 


40 


Ro(S d \A) — o-o xyz ^ ~ 1 ^ 


= 


h = 
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-continued 

2, if and € D n £S t d \A), n B € {n b h 

A set of base nodes in Gt(Si d \A) 

1, otherwise 


1, if € e {fCS/'M)} 

0, otherwise 


1, if the selected best pose of S t d \A {^*{S^\A)} 
0, otherwise. 


The estimate of relative assembly cost for S Z ^|A, 
R(S Z ^| A ), is the addition of the estimate of the relative 
assembly cost due to required reorientations, R 0 (S Z ^| A), 
and the estimate of the relative assembly cost due to 
required part matings, Rs(S z ^| A): 

Since the reorientations considered in the above compu- 
tation of R 0 (S Z ^|A) allow part matings to be done by 
stacking operations, R^S/^A) can be estimated simply 
as 


R£S t d \A)=/3o\-N, 

where N represents the number of parts included in 
S Z ^|A. 

An alternative way of estimating R(S Z ^|A) is that 
R 0 (S Z ^| A) is estimated based only on the required num- 
ber of base flip-overs, so as to avoid costly reorienta- 
tions as much as possible. Instead, the assembly of S z ^| A 
allows non-stacking operations the cost of which should 
be incorporated into R^S^I A). The computation of 
R 0 (S Z ^| A) due only to the base flip-overs can be com- 
puted simply by applying different a^ to the above equa- 
tion: 




I, ^ +£ and € D n £S d \A), n B € {n b }: 

A set of base nodes in G/£St d \A ) 

1, otherwise 


65 


where 


Whereas, R^S^I A) can be computed by 

Rs(Sf\A)=N.lJ3oi - U% mayk~ 

S?\A))l 

where 

tW | (U e (S/|A), £=x, y and 2 }. 

Then, R(S Z ^[ A) is chosen as the lower estimate from 
the two alternatives described above. 

Finally, the global estimation of the relative assembly 
cost, R, due to directionality, can be obtained by 
R=R(S^|A)+R(A-S^|A). 

5.3 Manipulability 

A subassembly subject to either a reorientation and- 
/or a translation for mating should be easily manipula- 
ble by tools or hands. The term manipulability of S z ^| A 
is used to quantify the efficiency in orienting S z ^| A and 
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in handling of S^| A. The manipulability of Sp\A is 
closely linked to the size, shape and weight of S^| A. 
More specifically, the orientation efficiency can be mea- 
sured based on the symmetry and marked polarity in the 
geometry and weight of S^[A, whereas the handling 5 
efficiency can be measured based on the regularity in 
the size, weight and shape of S/*| A, and the flexibility 
and fragility of Sp\A which determine the need for 
special tooling. 

The manipulabilities of S Z ^|A and A— Sp\A affect 10 
the relative assembly cost of the local mating operation 
between Sp\ A and A — S/| A, since they directly influ- 
ence the relative assembly cost for the required reorien- 
tations as well as the mating motion. 

To take this into consideration, in the previous sec- I 5 
tion, the relative assembly cost for a reorientation, a, as 
well as the relative cost for a mating motion, ft were 
determined by multiplying the manipulability coeffici- 
ent, 7, to their nominal values, ao and /Jo- 

20 

TABLE 2 

The Criteria for Measuring Manipulability of a Part or a 
Subassembly 

Relative Assembly Cost 

Orientation Efficiency 25 

Part tangles, nests or shingles 5 

Asymmetric paxt without marked 5 

polarities of weight or geometry 

Asymmetric part with marked 3 

polarities of weight or geometry 

Symmetric part 1 30 

Part delivered to the assembly station 1 

with a known orientation 
Handling Efficiency 

Large off center weight potentially 5 

causing loss of orientation 

Very large parts 5 35 

Very small parts 5 

Fragile 3 

Flexible 3 

Irregular shaped part requiring 3 

special tooling 

Easily handled part with standard 1 40 

(tooling cau handle more than 1 part) 

The manipulability coefficient of S z ^|A, y(S^|A), 
can be determined as the ratio between the sum of the 
score for each criterion of orientation and handling 45 
efficiency for S p | A and the sum of the nominal scores 
for each criterion of orientation and handling efficiency. 

6 Selection of Best Subassemblies based on AO* Algo- 
rithm 

As indicated earlier herein, the selection of Sp\A 50 
based solely on the relative assembly costs involved in 
the local mating operation between S/| A and A — 
Sp\A may not produce a globally optimal assembly 
plan. Therefore, we adopt the AO* algorithm with a 
properly defined evaluation function to search for a 55 
globally optimal or suboptimal plan. 

The search space to which the AO* algorithm is 
applied can be represented by an AND/OR tree. 

The decomposition of an assembly A in backward 
assembly planning implies the expansion of an AND 60 
node (representing an assembly A) into its OR children 
representing the alternative decompositions of A, 
{(Sp\ A, A — S^| A), i— 1, . . . , l}, and its AND grand- 
children {S/*|A and A— Sj**| A, for i = 1, . . . , 1} at- 
tached to indivisual OR children, as shown in FIG. 8. 65 
The AO* algorithm searches for an optimal solution 
tree by expanding those AND nodes of the current 
potential solution tree that are open to expansion, and 


by evaluating the next alternatives based on an evalua- 
tion function. 

An AND tree is an AND/OR tree every AND node 
of which has no more than one OR child. A potential 
solution tree is an AND tree having the minimum value 
for the evaluation function at the current stage of 
search, whereas a solution tree is an AND tree with 
leaves consisting of only single parts. 

The Local Cost, c?(nP), associated with an OR node 
np t represents the relative assembly cost incurred by the 
local mating operation between Sp\A and A— Sp\ A. 
c/n P) can be computed by the weighted sum of the 
following three components: 

1) The relative assembly cost due to the stabilization 
of Sp\A and A— Sp\A by using holding devices 
and/or reorientations, as described in Table 1. 

2) The relative assembly cost due to the reorienta- 
tions and translations required for mating between 
S^| A and A— S^|A, as described in Table 2 and 
Table 3. Note that this cost is linked to the direc- 
tionalities and best assembly poses for S Z ^|A and 
A — Sp | A, as well as the manipulabilities of Sp | A 
and A — S z ^| A. 

3) The relative cost of the special processes assigned 
to A, the parent node of n P. 

The Accumulated Cost, c^Tp) associated with an 
AND tree, Tp, represents the weighted sum of the 
following two components: 

1) The sum of c i(np) for all n P, n P T p. 

2) The depth of Tp defined by the maximum depth of 
n P for all n P, n P Tp where the depth of an OR 
node is measured in terms of the depth among OR 
nodes without considering AND nodes. 

The Local Heuristic Estimate, h^n, 0 ), associated with 
an OR node n p, represents an estimate of the optimal 
relative assembly cost to assemble S^] A and A— Sp\ A, 
and can be computed by the weighted sum of the fol- 
lowing components: 

1) The relative assembly cost, R, associated with the 
directional uniformity of Sp\A and A— Sp\ A, as 
defined in the previous section. 

2) The relative assembly cost, S, associated with the 
internal stability of Sp\ A and A— Sp\ A: 

S=d\x(S?\A)+x(A-Sp\A)] 

where x(S z ^| A) and x(A — S^| A) represent the internal 
stability of S z ^| A and A— S^|A, respectively, and are 
defined by 


XC Sp\A) - 


X (A - Sp\A) = 


the number of floating liaisons in Gi{Sp\A) 
the average degree of a node in Gi(Sp\A) 

the number of floating liaisons in Gi£A — Sp \ A) 
the average degree of a node in Gj_{A — Sp\A) 


and d represents the relative assembly cost due to a 
holding device. 

Then the evaluation function e/Tp) associated with 
an AND tree Tp simply becomes 

TP) - c a { TP) + 7 ih^Tp), 

where 7 ) adjusts the contribution of h£TP) to e/Tp) in 
relation with c a (T/). 

7 Assembly Planning with DFA Analysis 



5,442,563 


19 

As shown in the previous section, the evaluation of 
the local cost, c/(n,°), at an OR node, n P, is based on the 
detailed analysis of cfaP) in terms of the stability, the 
directionality, the assembly pose and the manipulability 
associated with the assembly of the children of n P, as 5 
well as the cost of special processes assigned to the 
parent of n P. 

The result of this analysis at each OR node of the 
search tree can directly be used for the identification of 
the assemblability of a product and for the evaluation of 10 
DFA criteria, which can be fed back to the designer for 
proper design evaluation and modification. The assem- 
bly planner developed here has both the capability of 
selecting an optimal assembly partial order as well as 
the capability of conducting DFA analysis, serving as a 15 
powerful tool for automating the DFA evaluation and 
modification cycle in concurrent engineering. DFA 
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implying that there is only one way to disassemble the 
product in the first step. Node 1 is expanded to it’s 
AND children, (2,3). Further expansion of Node 2 can 
be done easily because the grouping of B1 and B2 pro- 
vides only one alternative in selecting direct subassem- 
bly. At Node 3, the abstract liaison graph is shown in 
FIG. 10. In the drawing, a circle represents a node; a 
dash line represents a floating liaison; a thin solid line 
represents a firm liaison; a heavy solid line represents a 
rigid liaison; single hatching represents merging; and 
cross-hatching represents grouping. Note that S3 is a 
group node resulted from the special process forest. 

The algorithm to identify direct subassemblies is ap- 
plied to S2 and the result is shown in Table 4. There- 
fore, at Node 3, the system identifies two alternative 
direct subassemblies, C and S3. These two alternatives 
are represented as the OR 


TABLE 3 


A DFA Analysis Table for an OR Node, n/°, n,° ~ A — S/£|A} 


DFA Analysis Category 


DFA Criteria Details 


Stability 

Total 

The number of 

The relative 


relative cost 

holding devices re- 

cost due to the 


due to the 

quired for the sta- 

required hold- 


need to sta- 

bilization of S/*|A 

ing devices 


bilize Sp\ A 

and/or A — S/*| A 



and/or A — S/*|A 

The number of re- 

The relative 



orientations re- 

cost due to the 



quired for the sta- 

required reori- 



bilization of S/^j A 
and/or A — Sr| A 

enation 

Manipulability 

Total 

The manipulabil- 

y(S/*| A) and 

and Direc- 

relative cost 

ity factors 

y(A - S/IA) 

tionality 

involved in 

(Refer to Table 2 



mating be- 

for more details) 



tween S/|A 

The best assem- 



and A — S/^j A due to 

bly poses for S/*|A 



manipula- 

and/or A — S,^|A 



bility and 

The number of re- 

The relative 


directionality 

orientations 

cost due to the 



required for mat- 

required reori- 



ing between S,^J A 
and/or A — Sr j A 

entations 



The translati- 

The relative 



nal motion during 

cost 



mating 

due to the re- 



between Sp\A 

quired trans- 



and/or A — S/*|A 

latinal motion 
for mating 

Process 

Total rela- 

The list of special 

The relative 


tive cost for 

processes assigned 

cost of indi- 


the special 

to A 

vidual special 


processes 
assigned to 
A 

Total cost 
at n/° 


processes 


analysis performed during the process of computing the 
local cost, cjpaP) associated with an OR node, n P, n P — 
{S/*| A, A — Sp\ A}, is summarized into the DFA analy- 
sis table for n P, as illustrated in Table 3. 

Now, the analysis of DFA for a given product can be 
accomplished based on the DFA tables associated with 
all of the OR nodes of the solution tree. 

Example: The Flashlight Assembly 

The AO* algorithm with the cost and heuristic func- 
tions defined in the previous section is applied to the 
flashlight assembly for finding an optimal solution tree 
and performing DFA analysis. FIG. 9 illustrates first 
several nodes of the AND/OR search tree formed by 
the AO* algorithm, where DFA analysis tables are 
attached to individual OR nodes. At root node(Node 0), 
which corresponds to the whole product of flashlight, 
there is only one direct subassembly (Node 1 ). SI or S2 


TABLE 4 


Valid Cut-Sets or Direct Subassemblies generated from 
G //A) for S2 in the Flashlight Assembly 

Cut-set 

Pj A 

A - P| A 

Result 

{l9»ho) 

C 

S3+G 

Valid 


S3 + G 

C 

Valid 

{l&iio) 

G 

S3 + C 

Failed in the LM test 


S3 + C 

G 

Failed in the LM test 

{lS’b) 

S3 

C+G 

Valid 


C + G 

S3 

Valid 


nodes, 5 and 6, in FIG. 9, while the two OR nodes are 
expanded to their AND children, (9,10) and (11,12), 
respectively. The system then calculates the evaluation 
function at Node 5, e/Node 5), and the evaluation func- 
tion at Node 6, e/Node 6), based on the local costs, 


55 


60 


65 
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c^Node 5) and c^Node 6), and the local heuristic esti- 
mates, h e (Node 5) and h^Node 6), as follows: 

At Node 5, the system identifies that 

1) C is self-stable, but S3 + G needs one holding de- 
vice due to Is- 5 

2) The best assembly pose of S3+G needs one reori- 
entation so that BC + BL+R is on top of G. The 
assembly pose of the product is initially given in 
such a way that C is on the right side to BD. 

3) C needs one reorientation so that S3+G can be 10 
stacked onto C. 

4) S3+G + C need to be reoriented after the mating 
operation of S3 + G and C. 

5) The manipulability coefficients for S3+G(yi), 
C(y2) and S3 + G+C(y 3 ) are assigned to 0.8, 0.4 15 
and 0.4, respectively, based on Table 3. 

Therefore, 


c c (Node5) — 15(one holding device) 

+ 10 X 0.8 + 10 X 0.4 + 10 X 
0.4(3 reorientations) 

+ 1 X 0.8(stacking operation) 

+ 2(depth of Node 5) 

= 33.8 
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25 


For the calculation of the local heuristic estimate at 
Node 5, the system identifies the following: 

1) The estimates of the relative assembly cost due to 
the directional uniformity, R, and the internal sta- 30 
bility; S, of C are zero, since C is a single part. 

2) S34-C+G consists of 4 parts with the maximum 

uniform directionality of 1 in + x. Therefore, 
R 0 (S3-bC + G)= 10X 1 = 10 and 

RXS3-hC+G)= 1x4=4. As a result, the estimate 35 
of the relative assembly cost due to the directional 
uniformity, R, of S3 + C4-G is 14. 

3) S3 + C+G has 3 floating liaisons and has the aver- 
age degree of node of 2. Therefore, the estimate of 
the relative assembly cost due to the internal stabil- 40 
ity, S, of S3+C+G is 22.5 (S = 15x3/2=22.5) 
where 15 is used for the relative assembly cost for 

a holding device. 

4) The effect of Node 5 on assembly parallelism, w p , 
can be estimated as 6 (w^= 10x3/5 = 6 with 10 45 
assigned as a weight). 

As a result, we have that he(Node 5)=42.5. There- 
fore, we have that 


<?/Node 5) = c^/Node 5) + ij/z/Node 5) 

— 76.3(with 7 } — 1.0). 

At Node 6, the system follows the same steps that are 
used for calculating e/Node 5), and results are the fol- 55 
lowing: 

Co(Node 6) =22.6 


Therefore, we have that 

60 

e/Node 6) = c c (Node 6) + 7)^e(Node 6) 

= 54.6(with 7 } = 1.0). 

Finally, comparing e/Node 5) and e/Node 6), the 65 
system selects Node 6 for further expansion. The result 
of such an expansion is shown in FIG. 9. This process is 
continued until all nodes become single parts. 
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FIG. 11 illustrates a process implementing the fore- 
going, in which the merging process is performed 
(block 110 of FIG. 11 ), followed by the grouping pro- 
cess (block 120), both of which reduce the number of 
potential direct subassemblies which may be identified. 
All direct subassemblies are identified (block 130 ). 
Then, an AND -OR tree is constructed from a knowl- 
edge of all direct subassemblies (block 140 ). The local 
cost of each OR node is computed (block 150 ), and the 
tree is “pruned” to eliminate the higher cost OR nodes 
(block 160 ). Finally, an assembly sequence is carried out 
in accordance with the pruned AND -OR tree. 

FIG. 12 illustrates a system for carrying out the in- 
vention, including a data base representing the assembly 
design definitions 180 (such as the liaison graph or map, 
the assembly constraints and special process con- 
straints), a computer 190 which performs the process of 
FIG. 11 to generate a sequence of assembly instructions 
200. 

8 Conclusion 

This specification contributes to automatic assembly 
planning closer to reality by 

1) Developing a backward assembly planner which 
handles the case where an assembly sequence is not 
the same as the reverse of a disassembly sequence. 

2) Improving planning efficiency with the reduction 
of search space by merging and grouping parts 
based on interconnection feasibility and special 
process precedence constraints. 

3) Establishing assembly process planning by incor- 
porating such special processes providing as test- 
ing, cleaning, etc, in assembly planning. 

4) Providing subassembly evaluation criteria with a 
direct connection to assembly cost: especially the 
method for the stability and directionality of an 
assembly is presented through which the required 
number of holding devices and reorientations dur- 
ing assembly is identified. 

While the invention has been described in detail by 
specific reference to preferred embodiments thereof, it 
is understood that variations and modifications may be 
made without departing from the true spirit and scope 
of the invention. 

What is claimed is: 

1. A method of selecting and performing a sequence 
of sub-assembly operations to construct an assembly 
said assembly containing a plurality of assembly parts, 
said method comprising: 

identifying all direct subassemblies of said assembly, 
said identifying comprising: 

(a) identifying potential subassemblies of said as- 
sembly, 

(b) for each potential subassembly, determining 
whether said potential subassembly is a direct 
subassembly of said assembly by: 

(1) determining whether there is a common axis 
of separation between said potential subassem- 
bly and the remainder of said assembly, 

(2) determining whether said potential subassem- 
bly is stable, and 

(3) determining whether there is a direct unob- 
structed path between said potential subassem- 
bly and the remainder of said assembly; 

representing said assembly as a tree of AND and OR 
nodes beginning with all direct subassemblies as 
bottom nodes of said tree and concluding with said 
assembly as the top node of said tree, said AND 
nodes corresponding to required decompositions of 
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said assembly, said OR nodes corresponding to 
alternative decompositions of said assembly; 

computing a local cost of each OR node of said tree, 
said computing comprising at least one of the fol- 
lowing steps: 

(a) determining the number of reorientations of said 
assembly parts required to construct said direct 
subassembly given initial orientations of said 
assembly parts, and 

(b) determining the number of mating operations 
required to construct said direct subassembly; 

computing a heuristic cost of each OR node corre- 
sponding to said local cost comprising computing a 
cost of holding said direct subassembly multiplied 
by a cost of internal stability of said direct subas- 15 
sembly; 

combining said local cost and heuristic cost of each 
OR node of said tree to form a combined cost; 

eliminating OR nodes of higher combined costs; and 

performing an assembly operation of subassemblies 20 
corresponding to said tree after the elimination of 
OR nodes of higher combined costs. 

2. The method of claim 1 wherein said identifying 
potential subassemblies comprises: 

identifying all liaisons connecting pairs of nodes of 25 
each potential subassembly of said assembly; 

for each of said liaisons, merging liaisons for which 
there is no accessibility for interconnection be- 
tween said subassembly and said assembly and 
merging liaisons which have no force-deliverable 30 
paths to said assembly. 

3. The method of claim 2 further comprising identify- 
ing special processes of said assembly of the type com- 
prising testing and cleaning which require certain nodes 
of said assembly to be processed together during a spe- 35 
cial process, and wherein said identifying potential sub- 
assemblies further comprises: 

grouping together nodes of said assembly which must 
be processed together during special processes of 
the type comprising cleaning and testing of subas- 40 
semblies, whereby to prevent identification of po- 
tential subassemblies comprising corresponding 
nodes of said grouping as direct subassemblies. 

4. The method of claim 3 further comprising identify- 
ing precedent constraints corresponding to required 45 
sequences of assembly of parts of said assembly, and 
wherein said identifying potential subassemblies further 
comprises: 

grouping together nodes of said assembly which must 
be assembled in a hierarchical sequence corre- 50 
sponding to said precedent constraints, whereby to 
prevent identification of potential subassemblies 
comprising corresponding nodes of said hierarchi- 
cal sequence as direct or subassemblies. 

5. The method of claim 4 wherein said determining 55 
whether the potential subassembly is stable comprises 
determining one of: 

(a) said potential subassembly has only one transla- 
tional degree of freedom and only one rotational 
degree of freedom, said one rotational degree of 60 
freedom being about said translational axis only, 
and 

(b) any floating or disconnected liaison of said poten- 

tial subassembly includes a node of the remainder 
of said assembly. 65 

6. The method of claim 5 wherein said determining 
whether said potential subassembly is stable comprises 
determining said potential subassembly is unstable if 
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either said potential subassembly is unstable or the re- 
mainder of said assembly is unstable. 

7. The method of claim 1 wherein said internal stabil- 
ity cost comprises a ratio of the number of floating 

5 liaisons in said subassembly and the average number of 
liaisons in said subassembly. 

8. A method of selecting and performing a sequence 
of sub-assembly operations to construct an assembly 
said assembly containing a plurality of assembly parts, 

10 said method comprising: 

identifying all direct subassemblies of said assembly, 
said identifying comprising: 

(a) identifying potential subassemblies of said as- 
sembly, 

(b) for each potential subassembly, determining 
whether said potential subassembly is a direct 
subassembly of said assembly by: 

(1) determining whether there is a common axis 
of separation between said potential subassem- 
bly and the remainder of said assembly, 

(2) determining whether said potential subassem- 
bly is stable, and 

(3) determining whether there is a direct unob- 
structed path between said potential subassem- 
bly and the remainder of said assembly; 

representing said assembly as a tree of AND and OR 
nodes beginning with all direct subassemblies as 
bottom nodes of said tree and concluding with said 
assembly as the top node of said tree, said AND 
nodes corresponding to required decompositions of 
said assembly, said OR nodes corresponding to 
alternative decompositions of said assembly; 
computing a local cost of each OR node of said tree, 
said computing comprising at least: one of the fol- 
lowing steps: 

(a) determining the number of reorientations of said 
assembly parts required to construct said direct 
subassembly given initial orientations of said 
assembly parts, 

(b) determining the number of mating operations 
required to construct said direct subassembly, 
and, 

(c) determining a relative cost of holding devices to 
stabilize said direct subassembly and the remain- 
der of said assembly, and determining the rela- 
tive cost due to reorientation of said assembly 
parts; 

computing a heuristic cost of each OR node corre- 
sponding to said local cost comprising computing a 
cost of holding said direct subassembly multiplied 
by a cost of internal stability of said direct subas- 
sembly; 

displaying to an operator a table for each of said OR 
nodes summarizing all of factors of said local cost 
and heuristic cost; 

combining said local cost and heuristic cost of each 
OR node of said tree to form a combined cost; 
eliminating OR nodes of higher combined costs; and 
performing an assembly operation of subassemblies 
corresponding to said tree after the elimination of 
OR nodes of higher combined costs. 

9. The method of claim 8 wherein said identifying 
potential subassemblies comprises: 

identifying all liaisons connecting pairs of nodes of 
each potential subassembly of said assembly; 
for each of said liaisons, merging liaisons for which 
there is no accessibility for interconnection be- 
tween said subassembly and said assembly and 



5 , 442,563 


25 

merging liaisons which have no force deliverable 
paths to said assembly. 

10. The method of claim 9 further comprising identi- 
fying special processes of said assembly of the type 
comprising testing and cleaning which require certain 5 
nodes of said assembly to be processed together during 

a special process, and wherein said identifying potential 
subassemblies further comprises: 

grouping together nodes of said assembly which must 
be processed together during special processes of 10 
the type comprising cleaning and testing of subas- 
semblies, whereby to prevent identification of po- 
tential subassemblies comprising corresponding 
nodes of said grouping as direct subassemblies. 

11. The method of claim 10 further comprising identi- 15 
fying precedent constraints corresponding to required 
sequences of assembly of parts of said assembly, and 
wherein said identifying potential subassemblies further 
comprises: 

grouping together nodes of said assembly which must 20 
be assembled in a hierarchical sequence corre- 
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sponding to said precedent constraints, whereby to 
prevent identification of potential subassemblies 
comprising corresponding nodes of said hierarchi- 
cal sequence as direct subassemblies. 

12. The method of claim 11 wherein said determining 
whether the potential subassembly is stable comprises 
determining one of: 

(a) said potential subassembly has only one transla- 
tional degree of freedom and only one rotational 
degree of freedom, said one rotational degree of 
freedom being about said translational axis only, 
and 

(b) any floating or disconnected liaison of said poten- 
tial subassembly includes a node of the remainder 
of said assembly. 

13. The method of claim 12 wherein said determining 
whether said potential subassembly is stable comprises 
determining said potential subassembly as unstable if 
either said potential subassembly is unstable or the re- 
mainder of said assembly is unstable. 
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